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I. INTRODUCTION AND OVERVIEW 

Understanding what differentiates the collective behavior of a large quantum system from that of its classical 
counterpart is a challenging problem of theoretical physics. An important example here is a phase transition, due to 
which the system becomes ordered in one or another way. Intuitively, it is clear that intrinsic randomness of quantum 
systems should enhance thermal fluctuations in their work against ordering. For the most of realistic models, the 
theoretical explanation of such quantum effects is typically based on additional simplifications and uncontrolled 
approximations. However, as we show in this article, for a certain class of important realistic quantum models 
a rigorous mathematical theory based on the so called 'first principles' can be developed in detail, qualitatively 
agreeing with relevant experimental data. These models are quantum anharmonic crystals, used, e.g., in the theory of 
structural phase transitions in solids triggered by the ordering of light particles. At low temperatures, the mentioned 
quantum effects prevent such light particles from ordering and thus stabilize the crystal. We call this effect quantum 
stabilization, see (). Its mathematical description is based on the use of infinite dimensional (path) integrals and the 
corresponding path measures, which serve here as mathematical models of thermodynamic phases. This allows us to 
apply powerful tools of modern mathematical analysis, which, however, makes the theory quite involved. The purpose 
of the present article is to explain the main aspects of this theory in the way accessible also to non-mathematicians. 
We do this by interpreting the corresponding ideas and results published in our recent monograph (). As we aim 
at reviewing the results, but not the literature, the quotations are restricted to a necessary minimum. The reader 
interested in the bibliography on this topic is cordially invited to find it in the Bibliographic Notes to each chapter of 
()• 

An anharmonic oscillator is a mathematical model of a point particle moving in a potential field with multiple 
minima and sufficient growth at infinity, which makes the particle motion confined to the vicinity of the point of 
its (possibly unstable) equilibrium. In the simplest case, the displacement q of the particle from this point is one- 
dimensional and the corresponding potential energy V(q) is a symmetric continuous function with two minima at 
±g», separated by a barrier. In this case, one speaks of double-well potentials, an example of which is given in Fig. 1. 
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Let now an infinite system of such point particles be arranged into a crystal. That is, each particle is localized 
in the vicinity of its own crystal site. Suppose also that the particles interact with each other. The corresponding 
model, called an anharmonic crystal, is often used in solid state physics to describe, e.g., ionic crystals containing 
localized light particles oscillating in the field created by heavy ionic complexes. The ordering of such light particles 
may trigger a structural change of the whole crystal, see (). A particular example of this sort can be a KDP-type 
ferroelectric with hydrogen bounds, in which such light particles are protons or deuterons performing one-dimensional 
oscillations along the bounds, see (). Another relevant physical object is a system of apex oxygen ions in YBaCuO-type 
high-temperature superconductors (). Anharmonic oscillators are also used in models describing the interaction of 
vibrating particles with a radiation field (), strong electron-electron correlations caused by the interaction of electrons 
with vibrating ions (), or charge transfer along hydrogen bounds (). 

If the particle motion obeys the laws of classical mechanics, the low energy states of the particle are degenerate. 
By this we mean the existence of different states with the same energy, as the particle is confined to one of the 
wells if its energy is less than the height of the potential barrier separating the wells. Then the low temperature 
equilibrium thermodynamic states (phases) of the corresponding anharmonic crystal can also be multiple. If this is 
the case, as the temperature changes the crystal undergoes a phase transition - a collective phenomenon caused by 
the interaction and by the degeneracy mentioned above. If the particle mass is sufficiently small, it moves according to 
the laws of quantum mechanics, and hence can tunnel through the potential barrier. This tunneling motion eliminates 
the degeneracy, which might affect the ability of the phase transition to occur or even suppress it completely. On 
the theoretical level, such 'quantum effects' were first discussed in (). Later on, a number of publications on this 
item appeared, see () and the literature quoted therein. The key conclusion of those works is that the quantum 
effects become strong whenever the particle mass m p h gets small. This agrees with the experimental data on the 
isotopic effect in hydrogen bound ferroelectrics () and in YBaCuO-type high-temperature superconductors (). At the 
same time, experiments show that high hydrostatic pressure applied to a ferroelectric crystal diminishes the phase 
transition temperature, see, e.g., () and the table on page 11 in (). Hence, the reduction of the distance between the 
wells amplifies quantum effects. From this one concludes that also the shape of the localizing potential V plays an 
important role here. 

In the mathematical theory of classical anharmonic crystals, a thermodynamic phase appears as a Gibbs measure 
defined by a specific condition, involving the potential energy of a single particle and the interaction between the 
particles. Usually, this condition is formulated as the Dobrushin-Lanford-Ruelle (DLR) equation which the Gibbs 
measures in question have to solve. For quantum infinite particle systems, thermodynamic states are usually described 
as linear functionals on non-commutative algebras of observables, see (). However, for quantum crystals the direct 
construction of such states appears to be impossible as the relevant observables are unbounded operators. In view of 
this, the results mentioned above were based on rather indirect indications, e.g., the lack of a symmetry breaking () 
or of an order parameter (). It soon had become clear that a mathematical theory, which describes phase transitions 
and quantum effects in a unified way, and which as well takes into account the role of the localizing potential, ought 
to be based on a precise definition of thermodynamic phases. Such a theory was elaborated in a series of works, see () 
and the citations therein. The key aspect of this theory is the use of path integral techniques in which thermodynamic 
phases are constructed as probability measures on path spaces (path measures) that solve the corresponding DLR 
equations. Thereafter, necessary and sufficient conditions are derived for a phase transition to occur expressed in 
terms of the particle mass m p h, the interaction strength Jo, and a parameter characterizing the potential energy V. 
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In the simplest case where 

V(q) = -aq 2 + bq 4 , and q* = y/a/2b, (I) 

cf. Fig.l, both conditions can be expressed through one and the same parameter Q — (4m p ix/9H 2 )qf Jq, and have the 
following surprisingly simple form: (a) Q > d6(d) (sufficient); (b) Q > h (necessary). Here 0(d) is a certain function 
of the lattice dimension d, such that d9(d) > 1 and dO(d) — > I as d — > +oo, see ([751) and ([74]) below. The parameter 
x < I is determined by the spectral properties of the single-particle Hamiltonian ([5]) with h = 0, see ([M|) below. 

A self-contained presentation of path integral methods in the statistical mechanics of systems of interacting quantum 
anharmonic oscillators is given in the monograph (), addressed to both communities - mathematicians and physicists. 
The reader can find here a collection of facts, concepts, and tools relevant for the application of functional-analytic 
and measure-theoretic methods to problems of quantum statistical mechanics. This includes, in particular, a complete 
description of the theory mentioned above, as well as of its far-going extensions and implications. In this article based 
on that monograph, we present and explain some of its details. 



II. THE MODEL AND ITS THERMODYNAMICS 
A. The model 

To be more concrete, we assume that the considered model describes an ionic crystal and thus adopt the ferroelectric 
terminology. This means that we study a system of light particles oscillating in the field created by ionic complexes, 
which we suppose to be fixed (adiabatic approximation). A substance modeled in this way might be a KDP type 
ferroelectric (), in which such particles are protons - ions of hydrogen, constituting hydrogen bounds. Each particle 
carries an electric charge; hence, its displacement from equilibrium produces a dipole moment, proportional to the 
displacement. Therefore, the main contribution to the interaction between two particles is proportional to the product 
of their displacements. According to these arguments, our model is described by the Hamiltonian 

x,y x 

Here the sums run through a crystalline lattice, which we assume to be a d-dimcnsional 'hypercubic' lattice Z d 
consisting of vectors x — (x%, . . . ,Xd), Xi € Z, and Z being the set of all (positive and non-positive) integers. The 
displacement q x of the particle localized near such x = (xi, . . . ,Xd) is supposed to be a ^-dimensional vector, i.e., 

q x = (?x , ■ ■ -,qx^) and 

V 

(qx,q v )=Y,^ ) 1y ) ( 3 ) 

3=1 

stands for the corresponding scalar product. Regarding the interaction intensities J xyi we suppose that there exists 
a real- valued function ip such that 

J xy = <p(\xi-yi\,...,\xd-yd\). (4) 

We say that the interaction has finite range if ip vanishes whenever either of its arguments exceeds a certain value. In 
particular, <p can be nonzero only if one of its arguments is equal to one and all other arguments are equal to zero. 
This is the case of a nearest neighbor interaction. To avoid unnecessary mathematical complications, in this article 
we suppose that J xy has finite range. Note, however, that in () we consider the general case with the only assumption 
that 

:= ^2 \Jxy\ < 00, (5) 

see also (). In addition, we suppose that J xy > 0, which means that we consider ferroelectric interactions only. The 
single-particle Hamiltonian in ([2]) has the following form 

H x = ^-\ Px \ 2 + V(q x )-(h,q x ) (6) 



= H^ + U(q x )-(h,q x ), 
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where h <G M." is an external (electric) field, 



is the Hamiltonian of a harmonic oscillator of rigidity a > 0, and 

E%r) = V(&)-||&| 2 , (8) 

is the anharmonic part of the potential energy V. For the sake of convenience, we exclude from the latter the external 
field term. One observes that H x as given in ([6]) is independent of the choice of x; thus, the model (|2|) is invariant 
under the translations of the lattice T, d . The mass parameter m in (O and (JTJ) includes Planck's constant, i.e., 

m = TOph/fi 2 . (9) 

The Hamiltonian ([6]) is defined as a linear operator acting in the 'physical' complex Hilbert space H x = i 2 (K ly ) of 
square-integrable 'wave' functions, which have the standard physical interpretation. In particular, for ij) € % x such 
that 

U\\u x --= I \^{u)\Hu=l, 
and for a bounded linear operator A : T-L x — > 'H x , the (complex) number 

co^A) := {iP,Ai>) Hx = f ${u){A^){u)du (10) 

JR" 

is said to be the value of A in vector state uty. Here 

{<t>,i ) )u*= I <j>{u)ip{u)du, <j),ipeH x , 



is the scalar product in T-L x , and <p stands for the complex conjugate. In general, a state, w, is a map defined on the 
algebra £ x of bounded linear operators A : T-L x —> % x with values in the set of complex numbers C such that, for all 
<Ti,<T2 € C and all operators A\,A% e <L X , 

uj{<;iAi + <; 2 A 2 ) = ftw(Ai) + S2u{A 2 ), (11) 

uj(I) = 1, u(A) > for A>0. 

Here / is the identity operator and A > means that (i/>, Aip)-u x > for all V 1 € %c< Let wi and 012 be two states. 
For <j G [0, 1], we set 

w = + (1 - ?)w 2 . (12) 

Then uj is also a states, which is a convex combination {mixture) of wi and W2- The combination in (|12[) is called 
nontrivial, if a; 7^ uj\ and ui ^ u>2- By definition, a pure state is a state which cannot be decomposed into a nontrivial 
convex combination of other states. Vector states ([TO]) are pure, see Theorem 1.1.15, page 26 in (). 
In the sequel, we mostly consider the following versions of the model ©, ©: 

(i) the dimension v is arbitrary and V has the form 
cf. ©, 

V(q) = -a\q\ 2 + b\q\\ a, b > 0; (13) 

(ii) v = 1 and V is an even continuous function 
such that, for \q\ < \q\, 



W(q) - V(q) < W(q) - V(q), 



(14) 
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where 

r 

W(q) = ~b iq 2 + Y,biq 2 \ (15) 



1=2 

6iA>0, k > 0, for / = 2,...,r- 1, 



and r > 2 is a finite integer; 
(hi) ^ = f and is continuous, asymmetric, i.e., 
V(q) 7^ V(— q), and such that the following 
estimate holds 



V{q) > -C + D\q\ 2+e , qeR, (16) 



where C, D, and e are positive constants. 
Assuming that V in (|14j) is diffcrcntiable, we can rewrite this condition in the form 



-&i + J2 Ih^' 1 > for all > 0, 



1=2 

•2\ 



where v is such that V(q) = v(q ). Along with the function W we also use 

(20 



$ w = E 2 i-i ( z-i). 6,tfl - * - ' (17) 



;=2 



where r and all &z are the same as in (TTSl) . One observes that in cases (i) and (ii) with h — 0, the model is 0(f)- 
symmetric, that is, it is symmetric with respect to all orthogonal transformations of W, which constitute the groups 
0(v). Then the phase transition in the corresponding model can be associated with the breakdown of this symmetry. 
Case (hi) is essentially different from this point of view. 

We conclude this subsection by recalling that the Hamiltonian in corresponds to a quantum model. In particular, 
this means that p x and q x are (unbounded) operators, which in view of obey the following commutation relations 



P. 



U) q U') _ q p p U) = _ W j , i = V=l. (18) 



B. A single anharmonic oscillator 

As we shall see below, important thermodynamic properties of the model ([2]) are closely related to the quantum- 
mechanical properties of the anharmonic oscillator described by the Hamiltonian ([6]), in which we set h — 0. To 
avoid nonessential technical problems, we suppose that V in © is an 0(^)-symmetric polynomial and thus is of the 
following form 

r 

V(q) =^>M 2 ', c r >0, (19) 
1=1 

with r > 1, cf. (p~3]> . This, of course, includes also the case of the harmonic oscillator ([7]). For such V, there exists a 
subset T> of the physical Hilbert space H x such that H x with domain D is a self-adjoint operator, see Theorem 1.1.47, 
page 46 in (). Its spectrum consists entirely of real eigenvalues E n , n £ No, such that, for some E* > and n» € N, 
they obey the estimate 

E n > E„n 2r/{r+1 \ for all n > n*, 



1 Note that the only nontrivial transformation in O(l) is q i— > —q. 
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see Theorem 1.1.58, page 53 in (). In particular, the eigenvalues , n £ No, of H^ aT , are 

££ ar = (n + ^/2)A har , A har := V^T™, (20) 

see e.g. Proposition 1.1.37, page 41 in (). If v = 1, the eigenvalues -£„ are simple, which means that, for every n € No, 
the equation 

H x ip n = E n ip n 

has exactly one solution ij) n 6 The dynamical properties of the oscillator crucially depend on the following gap 
parameter 

A = inf (E n+1 - E n ) . (21) 

n£N 

Note that 

A ->• 0, as m ->• +oo. (22) 

The latter, in fact, is the classical limit h — > 0, see ([9]). It can be shown that, for a fixed value of the mass parameter 
m, 

lim (E n+1 - E n ) = +oo. (23) 

n— S-+00 

Thus, the infimum in (|21j) is strictly positive for finite values of m. As usual, for appropriate functions / and g, by 
writing / ~ g we mean that lim(/ j g) = 1 assuming that the meaning of the limit is clear from the context. In the 
harmonic case, the gap parameter is given in (1201) . According to Theorem 1.1.60, page 59 in (), A is a continuous 
function of the mass parameter m > such that, cf. (|22|) . 



A ~ A Q m- r / (r+1 \ as m -> 0, (24) 

where r is the same as in (|19|) and ([15]). The proof of these properties goes as follows. First, by the analytic perturbation 
theory for self-adjoint operators one shows that the eigenvalues E n , and hence the increments E n+ i —E n , continuously 
depend on m. This, however, does not yet mean the continuity of A, as the infimum of an infinite number of continuous 
functions need not be continuous. But, by (|23p we see that, for some no G N, 

E n +i — E n > Ek+\ — Ek 



for any n > no and any k < uq. This means that the infimum in (|2ip can be found by means of a finite number of 
increments, i.e., 

A = min (E n+1 - E n ) , 

n— 0,1, ...,no 

which yields the continuity in question. Finally, the asymptotics (|24|) is obtained by means of a rescaling and then by 
the asymptotic perturbation theory for self-adjoint operators. From this we conclude that the following parameter 

R m := mA 2 (25) 

is a continuous function of m and 

R n r,Alm-^- 1) / (r+1 \ as m -> 0. (26) 

Noteworthy, R m , J as given in ([5]), and the rigidity a as in (J7|), are measured in the same physical units, cf. ([9]). We 
also note that R m — > in the classical limit m — > +oo, sec (|93"1) below. In the harmonic case, i?„ r = a, see ([20)) . By 
analogy, we call R m the quantum rigidity of the oscillator described by H x . By the mentioned continuity and by (|26|) , 
for r > 2, R m takes any value from an interval [R^K +oo) for some R^ > 0. As we shall see below, for big enough 
values of R mi phase transitions in the model @ are suppressed at all temperatures. 
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C. Thermodynamic phases 

The Hamiltonian given in ([2]) has no direct mathematical meaning and serves as a formula for 'local' Hamiltonians 

H ^ = -\ 12 J ^(^' 9») + 12 Hx - ( 2? ) 

indexed by finite subsets A C 1 d . Let | A| stand for the number of lattice sites contained in A. Then the corresponding 
physical Hilbert space is Ha = -L 2 ((R")' A ' ). Among all linear operators acting in Ha, wc distinguish bounded operators 
A : Ha —} Ha- The set of all such operators Ca is a C*-algebra, see Sections 1.1 and 1.2 in () for more information 
on this topic. The local Hamiltonian defined in (|27p is an unbounded operator. However, by the assumptions made 
above Ha is a self-adjoint operator on some domain T>a C Ha- Moreover, one can show that, for any /3 > 0, the 
operator exp {—(3Ha) has finite trace, and hence the local partition function 

Za := trace exp (-(3H A ) , (28) 



is well-defined. Along with the Hamiltonian (|27|) we will use 'periodic' Hamiltonians defined as follows. Let Al, 
L e N, be the subset of the lattice 1 d consisting of all x = (xi, . . . , Xd) such that —L < Xi < L, for all i = 1, . . . , d. 
For Xi and ?/;, obeying the latter condition, we set 

\xi -Vi\h= min{|x l - Vi[,2L - \xi - y^}. 

Clearly, A^ is a cube, which one can turn into a torus by identifying its opposite walls. We denote it also by A^ and 
equip with the distance 



\x - v\l = \JW ~ Vi\l + --- + \x d - Vd\l- (29) 
After that, A^ becomes invariant with respect to the corresponding translations. Now we introduce 

J x V = <f{\ x i ~ Vi\l, ...,\xd- Vd\L), x, y e A L , 
where tp is as in (f5]). Then the periodic Hamiltonian in A^ is 

H L = -\ 12 Jxy(V*,%)+ 12 H *- ( 3 °) 
x,yeA L x£A L 

Like the Hamiltonian in (|2"T)) . it is a self-adjoint operator in the corresponding space Ha l , for which 

Zl := trace exp (— (3Hl) < oo. (31) 
The local Gibbs state qa at temperature T — l/ffks, ks being Boltzmann's constant, is defined to be 

qa(A) = trace {A exp (-/3H A )} /Z A) A e <£ A - (32) 
It is a positive linear functional on the algebra £a, of. (fTTj) . Moreover, it is the following mixture, cf. (|12[) . 

QA = -g- ^ eXp ("^A^) W <K 
n— 

of the vector states corresponding to the eigenvectors of H\ defined by 

In the same way, we define the periodic local state 

q l {A) = trace {Aexp{~pH L )}/Z L , A e £ A . (33) 

If not explicitly stated otherwise, all properties of the states qa that we mention in the sequel are attributed also to 
the periodic states ql. 



9 



By H0egh-Krohn's theorem, see page 72 in (), g A can be recovered from its values on the products of the form 

a£(Fi).-.(£ (*■„), nel, 

Fx, . . . , F n G $ A , ti, . . . ,i„ G M, 

where Sa is a 'rich enough' family of multiplication operators by bounded measurable functions F : IR' A ' — > C. Here 

af~[A) = ex.p(itH A )Aexp(-itH A ), A £ £ A , 

is the time automorphism : € A — > £a which in the Heisenberg approach describes the dynamics of the oscillators 
located in A. Thus, the Green functions 

G Fu ..., Fn (h,...,t n ) :=g A [a^F,) ■ ■ ■ ai(F n )] , (34) 

with all possible choices of n G N and then of JPi, . . . ,F n £ $ A determine the state g A . Each Green function admits 
an analytic continuation, also denoted by G?^ j? n) to the domain D n consisting of those (£i, ■ • • , Cn) G C", for which 

0<Im(Ci) < • • • < Im(Cn) < P- 

Furthermore, see Theorem 1.2.32, page 78 in (), Gf 1j ...,f„ is continuous on the closure D n and can uniquely be 
recovered from its restriction to the set 

:= {(Ci, . . . , Cn) G V n : Re(Ci) = • • • = Re(Cn) = 0}. 

Thus, the Matsubara functions 

rF 1 ,...,F„(ri,...,r n ) := G Flt ... >Fn (in,. . . ,iT n ), (35) 

with 

< n < r 2 • • • < r n < p, (36) 

and with all possible choices of n £ N and then of -Ft, . . . , F n G 3a i uniquely determine the state (J32J . By (|34|) . we 
can rewrite (|35p in the form 

r Fl ,..., F „(n, ...,r n ) = trace |> ie -(«-rO*Air aC -fr»-™)*A • ■ • F nC -(T»+i-r»)/r A j /Za> (37) 

where r n+ i = /3+rx and the arguments obey (1361) . The main ingredient of our technique is the following representation, 
see Theorem 1.4.5 in (), 

rir 1 ,...,F„(ri,...,r n ) (38) 

Fi(wa(ti)) ' ' ' F, l (wA(T„))^A(dwA), 

where z/a is the 'local Euclidean Gibbs measure' in A, which is a probability measure on the measurable space 
(n A ,B(n A )). Here 

Cl A := {lj a = (u x ) xeA : u x £ Cp}, 

C p := {0 G C([0,/3] -► R") : 0(0) = </>(J3)}, 

and C([0,/3] — > W) stands for the set of all continuous functions (paths) <f> : [0,(3] — > R", also called 'continuous 
temperature loops' in view of the periodicity property 0(0) = <f>(/3). Thus, Cp is the Banach space of paths. The 
space VL A is equipped with the corresponding product topology and with the Borel a-field B{VL A ). The measure v A 
has the following Feynman-Kac representation 

v A (du A ) = exp [-I A (u A )] XA{du A )/N A (h), (39) 
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in which 



1 f 13 
Ia(uk) = -~ J *y j (ux(r),uJy(T))dT 

-J2 (Koj x (r))dr + J2 U(oj x (T))dT 



(40) 



is the energy functional. We recall that U in (|4U|) is the anharmonic part of V, see ([5]), and that, for each r and 
x, lu x (t) is a ^-dimensional vector and thus (•, •) in (|40[) stands for the corresponding scalar product, cf. (|3]). The 
normalizing term N\(h) is defined by the condition 



N A (h) = / cxp [-Ja(wa)] Xa(^a), 



and the measure %a has the form 



(rfWA) = Yl X{dUx) 



(41) 



(42) 



Here \ is the H0egh-Krohn measure - a Gaussian probability measure on the Banach space Cp, constructed by means 
of the harmonic part of the Hamiltonian ([5]). It is defined by its Fourier transform 



J c ex P jf A ° ) ( T )^' ) (r)^ x (dw x ) = exp ( - \Y,J* £ 



S(T, T')A (j) (T)A^ (T')dTdT' 



where A is in Cp and S(r, r') is the propagator corresponding to a single scalar (y — 1) harmonic oscillator, which, in 
fact, is the corresponding Matsubara function (|35j) for n = 2 and Fi = F2 = q x . That is, cf. (|37|) . ([7]), and ([20l) . 



SYr, t') = —r— trace 



fe exp (-|r' - r|7^ ar ) g a exp (-(/3 - |r' - r|)^ ar ) 



(43) 



2^m^ [1 - exp (-/3A har )] , 



exp (-|r' - r|A har ) + exp (-(fi - |r' - r|)A har ; 
see pages 99 and 125 in (). In a similar way, one defines also periodic local Gibbs measures, cf. (j3"0")l . 

v L (duJA L ) = ex P ("Ai )] Xa l ( d UA L )/N L (h), 

with 



N L {h) 



exp {-I L {uj Al )] xa l (duA L ) , 



(44) 



(45) 



and 



E 

x,y£A L 



Jxy I (u} x (r),u} y {T))dr 





£ / (h,uj x (r))dr+ J2 / ' 



f/(w x (T))dr. 



Thus, the representation (|38|) leads to the description of the states (f32|) . <|33P in terms of Gibbs measures, similarly as 
in the case of classical anharmonic crystals. Here, however, each ^-dimensional vector q x is replaced by a continuous 
^-dimensional path lj x , which is an element of an infinite dimensional vector space. Going further in this direction, one 
can define 'global' Gibbs states of the model ((2J as the probability measures on the space of 'tempered configurations' 
satisfying the DLR equation, see Chapter 3 in (). It can be shown that the set of all such measures, which we denote 
by Q, is a nonempty compact simplex with a nonempty extreme boundary ex(^), the elements of which correspond 



11 



to the thermodynamic phases of our model, see also Chapter 7 in (). Thus, by a thermodynamic phase we understand 
an extreme element of Q. The latter means that such a measure cannot be decomposed into a nontrivial convex 
combination of other elements of Q. By virtue of the DLR equation, the set Q can contain either one or infinitely 
many elements. Correspondingly, the multiplicity (resp. the uniqueness) of thermodynamic phases existing at a given 
value of the temperature means that |ex((?)| > 1 (resp. \Q\ — 1). In general, there exist several ways of establishing 
whether the set of Gibbs measures is infinite or containing a single element. One of them consists in showing that 
there exists an element of Q 1 which is 'less symmetric' than the corresponding Hamiltonian, i.e., that a symmetry 
breaking occurs. Another way is to show that Q contains an element, which is not in cx(Q). Since the latter is a 
subset of Q, then Q is not a singleton, and hence \Q\ = +oo. The latter way will be discussed below in more detail. 



D. The free energy 

A more 'traditional' way of establishing phase transitions is based on the use of thermodynamic functions. In our 
case, this will be the free energy. Recall that we assume J xy to be of finite range. Along with the energy functional 
(1401 we consider also 

Ja(wa|0 = Ja(wa) (46) 



Jxy / (Wx(T),£ y (T))dT, 

V ..rlc JO 



l£A, y£A 

where Ja(wa) is as in (j40|) and all £ y with y £ A c := Z d \ A, belong to Cp. The last term in (|46| describes the 
interaction of the particles in A with some 'environment' defined by the configuration £ fixed outside A. We recall 
that I A (ui A ) depends on the external field parameter h £ R d . For xa as in (|4"21) . we set 

N A (h,£) = [ eatp[-7 A (w A |0]XA(dwA), 
Jn A 

which is well-defined as J xy has finite range, and hence the sum in (|46l) is finite. Obviously, if £ y = for all y E A c , 
then the above Na(H, £) coincides with N\(h) given in (|4T|) . We also note that the partition function (|28|) and N\(h) 
are related to each other by 

Z A = N K {h)Z\ a \ (47) 

where Z^ ar is the partition function of the system of non-interacting harmonic oscillators located in A. Explicit 
calculations yield, see page 139 in (), 

ZY = [traceexp(-/3^ ar )] |A| 



v\A\ 



1 _ e -0A*> 

where A har = yjajm is the same as in (l20l) and p3|) . In a similar way, we have 

Z L =N L (h)Z l r, (48) 

where Zl and Al(/i) are as in (|3"Tj) and (|45|). respectively. Therefore, NlQi), N\(h), and Aa(/i,£) are relative 
partition functions, and the latter one corresponds to the system of particles in A interacting with the 'environmental' 
configuration £. Let us now define 

F A (h,Z) = -^ x -UN A (h,Z), (49) 
F A (h) = -^-]nN A (h), 



and 



F L (h) = -p±-lnN L (h). (50) 
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Note that F A (h) = F A (h,0) and \Al\ = {2L) d . The functions just introduced are called local free energy densities. 
For a Gibbs measure \i £ G, we then set 

KW= f F A (h,ZM<%). (51) 
Ju 

This is the free energy density of the system of particles in A, interacting with the environment averaged with respect 
to fi. In the following, we are interested in the thermodynamic limits of the functions (j49|) - ([5"T]). To define such 
limits, we need one more notion which we introduce now. Let A be a finite subset of Z rf . We say that y £ Z d is a 
neighbor of A if: (a) y is not in A, i.e., y £ A c ; (b) there exists x £ A such that x and y are neighbors, i.e., 



\x-y\ = y/\xi - yi\ 2 H h \x d - y d 



2 



By OA we denote the set of all neighbors of A, |9A| standing for their number. A sequence {A„}„ e pf of finite subsets 
of 1 d is said to be a van Hove sequence, cf. page 193 in (), if: (a) A„ C A„+i for every n £ N; (b) for every y £ Z d , 
one finds n such that y £ A„; (c) lim n _>+oo |9A„|/|A„| = 0. An example of such a sequence can be the sequence of 
cubes Al„ for any strictly increasing sequence of positive integers L n . For a van Hove sequence {A„}„ 6 n, we denote 

F(h) = lim F An (h), F»(h) = lim Ft 1 (h), (52) 

where F An (h) and F£ (h) are as in (|49|) and (jSTj) . respectively. Furthermore, for a strictly increasing sequence {L n } 
and FL n {h) as in (fSTJ)) . we set 

F peI (h)= lim F in (/i). (53) 

n— f +oo 

It is known, see Theorem 5.1.3 on page 268 in (), that, for any [i£Q and for any van Hove sequence {A„} ne N, both 
limits in (|52p exist, do not depend on the choice of {A„}„ g N, and are equal to each other. Furthermore, for any 
sequence {L n } n eN such that L n — > +oo, the limit in ([53|) exists, is independent of the choice of {£ n }neNi an d is equal 
to the limits in ^2$. That is, 

F pcr (h) = F(h) = F"(h), for all fi£Q. (54) 

Therefore, F(h) is a universal thermodynamic function characterizing the model, that along with h depends also on 
/?. Note that in (), instead of F(h) we considered the pressure p(h) = —f3F(h). 

We recall that the model (JU with v > 1 is considered only with V as in (i). In this case the free energy density 
depends only on the norm of the external field, which we can choose therefore in the form (h, 0, . . . , 0) with h £ M.. 
In the remaining (ii) and (iii) cases, we have v = 1 and hence ftgl. Thus, from now on we assume that ft 6 1 in all 
cases. By (gTJ, and ([i7]l. we obtain 

- I^ Fa(/i) = ~w\ ^ Zh (55) 

Here M\(h) is the /ocaZ polarization averaged over A. Likewise, 

-s F 'C» = -OTS Zi (56) 

Here we have taken into account that the periodic state ql is invariant with respect to translations of the torus A^. 
The formulas just obtained relate the local free energy densities with the local polarizations. However, the convergence 
in (|52]l and (|53|) need not yield the convergence of the derivatives, i.e., of the local polarizations. Computing one more 
derivative in (JSHJ) we obtain, see page 220 in (), 



. x eA jQ 



lj x {t)) dr 



v A {dujA)vh{dLu A ) > 0. 



(57) 
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The latter inequality follows from the fact that we integrate a nonnegative expression, i.e., [■ • • ] . From this inequality 
we see that F\ is a concave function of ft. The same is true also for F^. Thus, the limiting free energy density given 
in (j52")l is also a concave function of ft. Concave functions always have one-sided derivatives 

Ai ± (ft) = -lim nh±e)-F(h) 

If M + (h) — M_(ft), then F(h) is differentiable at this ft, and one can speak of the global polarization 

M{h) := M+{h) = M_(/i) = -F'(h). (59) 

Given h* E K, we say that M(h*) does not exist if M + (ft*) ^ M_ (/i*). By the concavity of F, if the global polarization 
M(h) exists for a given h, it can be obtained as a limit of the sequences of local polarizations (|55| . (|56|) . which is 
independent of the choice of the corresponding van Hove sequences {A„} or {L n }. That is, 

M(h)= lim M A Jh)= lira M L (h). (60) 

Suppose now that, for a given ft* 6 R, -F(ft) is not differentiable at ft*, but is differentiable on the intervals (ft.* — e, ft*) 
and (ft*, ft* + e) for some e > 0. Then M + (ft*) ^= M_(ft*), i.e., the polarization M(h) is discontinuous (makes a 
jump) at ft*, which physicists interpret as a phase transition. However, we failed so far to prove that the discontinuity 
of M(h) implies the multiplicity of phases, i.e., |ex((?)| > 1, even for the simplest versions of the modeH (0). The 
converse statement holds true in the following form, see Theorem 5.3.3, page 276 in (). 

Theorem 1 For v = I, if the free energy density is differentiable in ft at a given ft* and (3, then the set of Gibbs 
measures Q is a singleton, and hence there exists only one thermodynamic phase at these (3 and ft* . 

In the next subsection, we return to the connection of the properties of F(h) with phase transitions. 

We conclude this subsection with the following remark. In some cases, it is possible to find values of ft where F(h) 
is differentiable for all f3. In particular, see Theorem 5.2.3 on page 273 in (), if v = 1, 2 and V is as in (|13| . then F(h) 
is infinitely differentiable at all ft 7^ 0, and hence M(h) exists at all such ft. This result follows from a generalization 
of the celebrated Lee- Yang theorem, see also Theorem 5.2.4 on page 275 in (). 



E. Phase transitions and order parameters 

We recall that according to our definition a phase transition in the model ((2|) occurs if it has multiple thermodynamic 
phases existing at the same values of /3 and ft. Mathematically, this means that the set of Gibbs measures Q, and 
hence its extreme boundary ex(C?), contain more than one element. So far, for models like ([2]) there exists only one way 
of establishing the latter fact. Namely, one shows the existence of an element of G, which is not extreme (not a pure 
state). Recall that ex(Q) is a subset of Q. A candidate to be such an element is the so called 'periodic state' which 
one obtains as the limit of the measures (|44|) as L — > +00. For ft = and big enough j3, this state can be nonergodic 
with respect to the translations of Z d , whereas extreme states are always ergodic. We use this argument in the next 
section. However, this way essentially employs the 0(^)-symmetry, and hence does not cover case (iii). Therefore, to 
mathematically describe phase transitions in a wider class of models one needs other definitions, consistent with the 
physical point of view on this subject. In this subsection, we present and analyze such alternative definitions. 

The first definition is based on the celebrated L. D. Landau classification, which employs the differentiability of the 
free energy density in ft, see, e.g., Chapter I in (). Of course, this property depends also on the value of /?. Namely, 
we say that the model (J2j) has a first-order phase transition at certain values of /? and ft* if at this f3, the polarization 
M{h) is discontinuous at ft*. The model has a second- order phase transition at (3 and ft* if the first derivative F'(h) 
is continuous, but the second derivative F"{h) is discontinuous at ft*. Note that these definitions can be applied to 
version (iii) of our model. 

The second definition employs an order parameter. It can be applied to cases (i) and (ii) of the model ([2]) with 
ft = 0, which we suppose to be 0(^)-symmetric. One can show that the local states (|3"2"j) and (|3"3")l are such that 



2 Such a statement holds true for the Ising model, see Theorem III. 3. 11, page 260 in (), as well as for the model of a nonideal gas in M. d , 
see (). 
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PAd&l 2 ) an d Ql{\<1x\ 2 ) exist and are finite, in spite of the fact that Ig^ | 2 is an unbounded operator. Then we can 
consider the Matsubara function, cf. (|37|) and (l43l) . 



r^(r,r') (61) 

— > trace 
7,t ^ 



q U) e -(r'-r)H L q U) e -(/9+ T -r')Hr 



Here < r < t' < (3, x and y are points in the cube A^, L £ N, and Hl and Zl are given in ([3"U|) and (|31|) . 
respectively. According to (pI5)) we have 

^y(r,r') =J2[ 4 j) (T)4»{r')Mdu> AL ), (62) 

j=l JU *L 

where vl is given in (|44[) . By the Cauchy-Schwarz inequality, and then by a property of the integrals with respect to 
vl , it follows that 

|r^(r,r')| < g L (\q x \ 2 ) < C. (63) 

Note that Ql{\ox\ 2 ) is independent of The meaning of the second inequality in (|6"3"|) is that the bound C is 
independent of A^. 

As a function of r and r', r^ y in (|62"|) can clearly be extended to all values of (r, r') € [0, /3] 2 . This extension is a 
continuous function of 

\t-t'\ p :=min{|r-T'|;/3- |r-r'|}. 
Then the following integral makes sense 



D L xy =p- 1 I ! T L xy {T,T')dTdT> . 

Jo Jo 



(64) 



If the model ((2]) is 0(^)-symmetric, as it is in cases (i) and (ii) with h = 0, D xy describes static correlations between 
the oscillators located at x and y. In view of the translation symmetry of vl, D xy is a function of the periodic distance 
\x — y\L, see (|29j). Now we introduce, cf. ([62]) . 

Ml Ll x,y£A L 

2 



If D xy decays sufficiently fast as \x — y\L — > +oo, then the double sum in the first line of (|65p is of order strictly less 
than | | 2 , and hence Pl — > as L — > +oo. On the other hand, by means of (|S3"|) we readily obtain that Pl < C, 
which means that the sequence {PlIlgn is bounded and hence contains convergent subsequences. The largest limit 
of such subsequences, i.e., 

P = lim sup Pl , (66) 

L— >-\-oa 

is called order parameter. We say that, at a given value of (3, there exists a long range order, if P > for this /3. In 
this case, we also say that the model is in an ordered state. Thus, we have the following three definitions of a phase 
transition: 

(a) employing thermodynamic phases: |ex(C?)| > 1; 

(b) employing the polarization jump: 

M + (h) ± M_(fc); 

(c) employing the order parameter: P > 0. 

The connections between these notions will be analyzed in the next section. 
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F. Infrared bounds 

The fact that D xy is a function of the periodic distance (f2T)f allows us to use the Fourier transformation in the 
torus Al. Let denote the Brillouin zone for A^, that is, the set consisting of p — (pi, ■ ■ -Pd) such that 

p,. = — 7r + nsi/L, Si = 1,2, ...,2L, i = l,,..,d. 

By (—Tr,ir] d we denote the set 

{p = (pi, . . . ,p d ) : p t € (-7T,7r]}, 
which is the Brillouin zone for the whole crystal. For p £ (— Tr,n] d , we define 

3 p = E D^exp(i(p,x-y)). (67) 

The latter sum is finite for finite L, but can be divergent in the limit L — > +oo. Traditionally, this is called the 
infrared divergence. The restriction of to h* L can be used to define the inverse transform 

D x V = T7— | Y] ®p exp(-i(>, x - y)) 

Suppose now that we are given a continuous function B : (— 7r,7r] d — > (0, +oo] with the following properties: 

(a) / B(p)dp < +oo; (68) 

J(-Tr,ir] d 

(b) D% < B(p), for all p^O. 

Since B is independent of L, by (b) one can control the infrared divergence mentioned above. That is why (b) is 
called an infrared bound, see () for more detail on this item. In view of (a), for x, y G Z d the following quantity 

B xy = I B(p) cos(p, x - y)dp (69) 



(-7r,7r] d 



is well-defined. For x, y e A^, we also set 



B x V = 7T—i Y\ B(p) cos(p, x-y), 



and B xy = if either of x,y belongs to K C L . One can prove, see Proposition 6.1.2 on page 284 in (), that, for every 
x and y, B x — > B xy as L — > +oo. Furthermore, see Lemma 6.1.3 ibid, for every A^ and any x, y <E A^, under the 
conditions in (l68l) the following holds 

D'xy > (D r xx - B xx ) + B' xy . (70) 

Note that both D xx and B xx are independent of x. The meaning of (|70|) is that the infrared bound (b) in (|68|) allows 
one to control from below the decay of the static correlations. Suppose now that there exists a positive $ such that, 
for any cube Al, 

D L XX > 0. (71) 
Suppose also that ([68| holds and that B xy obeys the following conditions: 

(a) d>B xx ; (72) 

(b) Km B xy = 0. 

\x—y\^-\-oc 
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If this is the case, for some e > 0, all and all x, y E Al such that \x — y\ > 5(e), by (|70)) we get that D^ y > e, 
which yields in (|65l) that Pl > s/j3 and hence P > 0, see (l66l) . Thus, (|71l) and (|72|) imply the existence of the long 
range order mentioned above. 

It turns out that, for our model (|2|) with any value of the dimension veN, and with any V, the function B(p) can 
be found explicitly. The only conditions are that the lattice dimension d should be at least 3 and the interaction in 
@ should be of nearest neighbor type, i.e., J xy = J > whenever |a; — y\ = 1, and J xy = otherwise. Then this 
function has the following form, see Corollary 6.2.9, page 303 in (), 

d 

£ (P) := - COS P*)' 
i=l 

Note that B(p) ~ C|p|~ 2 , as \p\ — > 0. Hence, it satisfies condition (a) in (|68| only for d > 3. By (|69l) . we get 

B xx = v6(d)/2J, (74) 

1 ' (2nr J { _ nMd e(p) 

Then, for any i} > 0, condition (a) in (|72p can be satisfied by taking big enough J. The validity of (b) in (|T2")) follows 
by the Riemann-Lebesgue lemma, see Proposition 6.2.10, page 303 in (). 



III. THE RESULTS 
A. Phase transitions 

In the first subsection below, we analyze the relationships between the three definitions of a phase transition given 
above, and then present the corresponding sufficient conditions for the particular versions of our model. 



1. The order parameter and the first-order phase transition 

First of all we note that the order parameter (|6"6")l can be defined without using any path integrals, see (|6ip . (|64l) . 
and (j6"5j) . The same is true for the free energy density (f53"|) since 

see (|5U)) . dH]), and (l3"Tj) . However, the path integral representation of these quantities allows us to apply here measure- 
theoretic tools, which proved to be useful in classical statistical physics. Namely, by means of Griffiths' theorem, see 
Theorem 6.1.7, page 286 in (), we obtain the following result. 

Theorem 2 Let the dimension v be arbitrary and the potential V be 0(y)- symmetric and such that the estimate H16\) 
holds. Then the one-sided derivative M+(0), see 158\) . and the order parameter P , see 166\) . obey the estimate 

M+(0) > VP. (75) 

Hence, if the model with such V is in an ordered state, then it undergoes a first- order phase transition. 

Indeed, by the assumed symmetry we have M + (0) = — M_(0). Hence, P > by (|75|) yields M + (0) > 0, and then 
M_(0) < 0. Note, however, that the positivity of the order parameter need not yet mean that the thermodynamic 
phases are multiple. To prove that this is the case we have to impose further restrictions on v and V . 
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2. Phase transition in the |g| 4 case 

Here we suppose that v is arbitrary and V is as in (fT3)l . Recall that a > is the rigidity parameter, see ©. Then, 
for a and b as in (fT3|). we set 



Note that #* > whenever a > and thereby the potential energy in ([6]) has multiple minima. For v = 1, V(q) has 
two minima at ±g*, with g* = v3$*> cf. ([I]). 
Let / : [0, +oo) — > (0, 1] be defined as 

/(0) = 1, and /(utanhu) = u^ 1 tanhu, for u > 0. 

Then, for a fixed 7 > 0, the function 

0(t):= *y/(t/7). ^>0 (77) 

is diffcrentiable and monotone increasing to 7 2 as i — > +00, that readily follows from the definition of /. Thus, we 
set 7 = 4mi?* and obtain that the function 

^-J<t>{P) = 2/3Jtf*/(/3/4rm?) 
Zm 

increases as f3 — > +00 and tends to 8mi? 2 J. Therefore, if the following condition holds 

8m^J > 9(d), (78) 

then the equation 

20J0*f(fi/4m&*) = 6>(d). (79) 

has a unique solution, which we denote by The following statement, see Theorem 6.3.6, page 308 in (), gives a 
sufficient condition for a phase transition to occur. 

Theorem 3 Let d > 3 7 the interaction intensity be J xy = J > whenever \x — y\ = 1, and J xy = otherwise, and 
let V be as in Suppose also that the condition J7ff| ) is satisfied, and hence J7ff| ) /ias i/ie solution /3*. Then, for 

P > P*, the following holds: (a) |ex(0)| > 1; (b) P > 0; (c) M + (0) > M_(0). 

The implication (6) (c) follows by Theorem [2j The proof of (a) and (b) relies on showing that the estimates (ffTj) 
and (|72p hold with as in (|73p and (I7H) . First, for the state (|33p and the Hamiltonian (f5D|) . and for an appropriate 
operator A, one readily gets that 

g L ([A, [H L ,A]])>0, 

where [Hl, A] = HlA— AHl. Applying this formula with A = p x \ x g and j = 1, . . . , v, and taking into account 
the form of V and the commutation relations (|18l) , we obtain 

QL(\q x \ 2 )>$*- (80) 
Next, we employ the Bruch-Falk inequality, see also () and page 392 in (), and obtain from the latter that 

d l > - (si) 



where / is as in (1771) . This estimate and ([71)1 yield 



For p > p m , the right-hand side of (|82|l is strictly positive. Then D% v does not decay to zero as L — > +00 and 
\x — v\l — > +00. This yields P > 0, sec (|65| . This also yields that the periodic state is not ergodic with respect to 
the group of translations of Z d , implying that Q contains elements which do not belong to ex(CJ), see Definition 3.1.26 
and Corollary 3.1.29, page 207 in (). By this we get that \Q\ > 1, which proves also (a). 
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3. Phase transition in the v = 1 symmetric case 

Here we assume that V(q) = V(—q) and ([T4")) is satisfied. In addition, we assume that there exists J > such that, 
for all x, y € 7L d obeying \x — y\ = 1, the following holds 

Jxy > J- (83) 

Let $ be as in (ITTl) . Then the equation 

- 2bi + $(i?)= 0, (84) 



has a unique solution > 0, since 6i > 0, see ([15)) . In this case, we have the following result, see Theorem 6.3.8, 
page 310 in (). 

Theorem 4 Let the model be as just described and let the condition in |7ff[ ) be satisfied with defined by {8$ 
and J obeying i83\) . Then, for every f3 > /3* defined by |7ff[ ), </ie following holds: (a) |ex(C?)| > 1; (b) P > 0; (c) 
M+(0) > M_(0). 

Again, as above (b) implies (c) by Theorem [3J Claims (a) and (6) are proven by comparing the considered model with 
the model described by the Hamiltonian 



x,y: \x~y\ — l x 



(85) 



where is as in (j!5p . Then we use the GKS inequalities, which hold in the considered case, see Theorem 2.2.2, page 
163 in (), and a comparison method developed in (). By means of these tools we prove that if (a) and (b) hold for the 
model as in (|85|) . they hold also for that of ([2]) with such V and v = 1. On the other hand, the properties of W allow 
us to prove that (|81[) holds for the model (|85"T) with defined by (j84[) . The infrared bound (|67|) with i? p as in ([73]) 
also holds since the interaction is of nearest-neighbor type. Then the proof of (a) and (6) for the model ([85]) follows 
as for Theorem [3] 



4. Phase transition in the v = 1 asymmetric case 

In case (iii) of the model ([2]), the only result concerning phase transitions which we managed to get so far is a 
statement that the model parameters and the inverse temperature /3 can be chosen in such a way that the polarization 
M(h) becomes discontinuous at certain /i*, i.e., the model has a first-order phase transition, see Theorem 3.4 of () 
and Section 6.3.3 of (). 

Theorem 5 Suppose that v = 1, d > 3, the interaction is of nearest-neighbor type, and V is continuous and such 
that M6\) holds, which includes the case ofV{q) ^ V(—q). Then, for every (3 > 0, there exist m* > and J* > such 
that, for all m > to* and J > </*, there exists h* € K, possibly dependent on m, /3, and J, such that the polarization 
M(h) becomes discontinuous at h = h*, i.e., the model has a first- order phase transition. 

Here we again exploit the fact, used in the proof of Theorem [3j that the existence of a nonergodic state implies 
= +oo, and hence a first-order phase transition. Recall that the existence of a nonergodic periodic state would 
follow from the fact that the static correlation function D^ y — [Mi(/i)] 2 does not decay to zero as L — > +oo and 
1^ ™ u\l +oo. The latter means that there exists a sequence of integer numbers {L„}„ e N such that L n — > +oo and 
the following holds 

D^-[M L Ah)] 2 >S, (86) 

for some S > 0, all x, y, and all n such that x, y € Al„. Note that in the symmetric case, we have Ml(0) = 0, and 
hence ([70]) can be used. Suppose now that there exists ho such that M(ho) = 0. Since M(ho) = limMj^/in), see ([601) . 
for such ho ([8B"|) would follow from the fact that 

D% > S. (87) 

Thus, if ([8"T| holds and M(h) takes both negative and positive values, then either M(ho) = 0, and thus ([B"6"[) yields a 
first-order phase transition, or M(h) is discontinuous, which yields the same. 
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To realize this scheme we crucially use the properties of our path integrals. As was mentioned above, the free energy 
density F(h) is a concave function of ft 6 1, see (|57)l . Then the set of such h where M + {h) ^ M_(/i), and hence 
M(h) does not exist, is at most countable. Let 1Z denote the set of all h £ M. for which M(h) exists. Then any interval 
(a, d) cK contains points of 1Z. A tedious analysis of the properties of the path integrals used in our constructions 
yields the following estimate of the free energy density (f5U|) 

-F L (h)>he + V(/3) + Z(m), (88) 

see Theorem 5.2.2, page 271 in (). This estimate holds for all h > 0, J > 0, and to > too, where too > is a certain 
quantity which may depend on j3. Here e > is a constant, and 'J and S are certain real-valued continuous functions 
of /3 and m, which can be calculated explicitly. As the right-hand side of (f88| is independent of L, this estimate holds 
also in the limit L — > +oo, for all h such that M(h) exists. By the mentioned concavity of F(h), for any such h > 0, 
we have 

M(h) > [F(0) - F(h)} /h 

>£+[F(0) + *(/3) + S(to)] /h. 

Since the right-hand side of the latter inequality tends to e > as h —> +oo, M(h) becomes positive for big enough 
h. Thus, there exists h+ > such that M(h) > for all h £ TZ D (h+, +oo). Likewise, one can show that there exists 
h- < such that M(h) < for all h £ 7Z n (— oo, h-). Next, some additional analysis of the path integral in the 
representation 

Ql(c£)= / K(0)] 2 ^ L (dw AL ), 

JflA L 

cf. UH), yields, see Lemma 6.3.9, page 312 in (), that, for every positive j3 and there exist positive mi and J*, 
which may depend on fi and but are independent of h, such that, for any cube A^ and any h £ K, and for all 
J > J* and to > toi, the following holds 

QL (q 2 x ) > (89) 

Then the proof of Theorem [5] follows along the next line of arguments. Fix any and take to > to* := max{rao; mi}. 
Then both (J88]) and |J89]) hold. Hence, in view of ([89]). we have that dSTJ> and ([82]) hold for such /3 and to, and for all 
L and /i. Then we increase J, if necessary, up to the value at which the right-hand side of (|82|) gets strictly positive. 
Afterwards, all the parameters, except for h, are kept fixed. By (f70|) . the positivity of the RHS of (|82|) yields (|87[). 
If M(h) were everywhere continuous, the fact that it takes both positive and negative values would imply that there 
exists ho such that M(/i ) = 0, which together with ([57]) would yield (|56")) . The latter, however, contradicts the 
everywhere continuity of M, see Theorem [T] 



B. Quantum effects 

1. The stability of quantum crystals 

To understand what makes a crystal stable or unstable let us first consider the scalar (y = 1) harmonic version of 
the model |2|), in which H x = flj ar , see (J6j) . For this model, the global Gibbs states, and hence the thermodynamic 
phases as mathematical objects, can be constructed only if the stability condition 

J < a = TO[A hal f 

is satisfied, see j5]) and (|20|) . In this case, \Q\ = 1 for all /3, see also Theorem |6] below . If Jo = a, the harmonic crystal 
becomes unstable with respect to spatial translations; for d > 3, the set Q is still a singleton, and Q — for d < 2, 
because of the divergence of integrals as in (|74|) for such values of d. In the anharmonic case, due to the assumption 
(fTB"]) global Gibbs states exist for all Jo > 0, and the instability of the crystal can be caused by the 'effective' change 
of the equilibrium positions of oscillators, i.e., by a structural phase transition. Thus, according to Theorem [2l a 
sufficient condition for a c?-dimensional quantum crystal, d > 3, to have multiple thermodynamic phases at a given /3 
is 

W*/ ( -j^r) > d6(d), (90) 
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which in the classical limit m —> +00, cf. (|9|), takes the form 

/3J tf* > d0{d). 

The latter condition can be satisfied by picking big enough /3. Therefore, the classical anharmonic crystal with d > 3 
always has a phase transition - no matter how small Jo is. For finite m, the left-hand side of (I90p is bounded by 
Am^l Jo, and the bound is achieved in the limit /3 — > +00. Thus, if 

Q := 4m^J < d9(d), (91) 

the condition (|9U)) will not be satisfied for any (3 > 0. Although (|9T)|) is only a sufficient condition, one might expect 
that the phase transition cannot occur at any /3 if the parameter Q is small. This effect could be called quantum 
stabilization. The following statement shows that this really occurs. 

Theorem 6 Let v = 1 and the potential V be as in U9\) . Then the set Q of the Gibbs states of the corresponding 
model is a singleton at all (3 > under the following condition 

Jo<R m - (92) 
Thus, in this case there is no phase transitions at all temperatures. 

This statement is an adaptation of Theorem 7.3.1, page 346 in (). Thus, (|9"2")l is the condition under which the 
mentioned quantum stabilization holds. In view of (126[) . it can be satisfied by picking small enough m. 
To relate (|92"T) with (|9Tj) . we use the following result, see Theorem 7.1.1, page 338 in (). 

Theorem 7 Let v and V be as in Theorem^ Then the rigidity parameter obeys the estimate 

* (93) 

This estimate clearly shows that R m describes the tunneling motion between the wells. For if z9* is big, cf. (1751 . then 
R m gets small. On the other hand, for small the bound (J93J) becomes inessential, and hence R m can be arbitrary. 
For such we lose the control of R m . This includes the case of convex V where i9» = 0, as it is for harmonic crystals. 
Thus, if (|92l) holds, then 



Q < x := Am$lR m < 1. (94) 
It is known, cf. Proposition 6.3.5, page 308 in (), that d6(d) > 1, and 

' < °( d ) < T~~T' (95) 



d- 1/2 w d-V 

for d > 4. Therefore, for the phase transition to be suppressed at all temperatures, it suffices that (|94| holds, whereas 
(|9"H is the corresponding necessary condition. By (|95"T) . the right-hand side of (f9"Tj) is close to one for big d. Thus, for 
such d the gap between these two conditions is close to 1 — x. 



2. Normality of fluctuations in the v > 1 case 

Unfortunately, we cannot extend Theorem [5] to the vector case v > 1. The main reason for this is that the 
eigenvalues of the corresponding H x are no longer simple. The only result of this kind which so far we can get is that 
the fluctuations of the displacements q x remain normal at all temperatures under a condition similar to (|92[) . 

Here we again suppose that the potential V is 0(i/)-symmetric. Then the operator which describes fluctuations of 
the displacements of oscillators in a given A is 

1 1 xeA 

Let us introduce the following parameter 

P A ] = \Y. f [^ F M F ^(T,T')dTdT>. 

P^Jo Jo J- A ^ 
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Then P L = P)^, cf. (65]). Clearly, if P > 0, then Pj^ -> +oo as L -> +oo, for any a G [1,2). We say that the 
fluctuations of the displacements of oscillators are normal if the Matsubara functions ([37]) for the operators (|96|) with 
fj = l remain bounded as A — > Z d . 

Since V is 0(V)-symmetric, there exists i > 0, such that V^g) = u(|g| 2 ). Let A be the gap parameter (|2"Tj) of 
the Hamiltonian 

H x :=^p 2 x + v(q 2 x ), (97) 

of a one-dimensional oscillator with the same m as the considered z/-dimensional oscillator. The following result is an 
adaptation of Theorem 7.3.5, page 350 in (). 

Theorem 8 Let v > 1 be arbitrary, V obey \lb}) , and let the condition A92\) be satisfied in which A is the gap parameter 
i21\) of the Hamiltonian |ff7[ ). Then the fluctuations of the displacements of oscillators in the corresponding model 
remain normal at all temperatures. 

The proof of this result relies on the so called scalar domination estimates obtained in Theorem 2.5.2, page 175 of (), 
see also (). 



3. Decay of correlations 

For simplicity, we suppose here that the model ([2j is of type (i), cf. (fl3| . with v = 1. For this model, the 
Matsubara function (|62|) describes the displacement-displacement correlations in the local periodic state (p3| . In the 
limit L — s- +oo, we obtain periodic Gibbs states, possibly multiple. Under the stability condition (|92|). there is only 
one such state. Then we set 

T xv (t,t')= lim r^(r,r'). (98) 

We know that T xy (r, r') should vanish in the limit \x— y\ — > oo, otherwise the corresponding state would be nonergodic. 
The following statement, which is an adaptation of Theorem 7.2.2, page 344 in (), describes the spatial decay of (|98|) . 



Theorem 9 Let the stability condition i9ty) be satisfied. Then the correlation function t98\) has the bound 

l xy [r,r)< to(A 2 + k 2 ) + T(p) dP ' (9Jj 



where 



T (p) = 2 ^2 Jx v \- sin (p> x -y)f 

y 



and JC is the set of all k = (2tt //3)k, k£Z. 

One observes that the right-hand side of (|99[) is the correlation function of a quantum harmonic crystal with the 
single-oscillator rigidity mA 2 . Thus, under the condition (|92j) the decay of (|98|) is at least as strong as it is in the 
corresponding stable quantum harmonic crystal. This result is valid for arbitrary J xy obeying ([S]). If the interaction 
has finite range, then T(p) ~ C\p\ 2 for some C > 0. In this case, the spatial decay of the right-hand side of and 
hence of T xv (t, r'), is exponential, see Theorem 7.2.4, page 345 in (). 



IV. CONCLUDING REMARKS 
A. To Section [TT] 

1. The model 

In ionic crystals, the ions usually form massive complexes which determine the physical properties of the crystal, 
including the instability with respect to structural phase transitions. Such massive complexes can be considered as 
classical particles which obey the rules of classical statistical mechanics. At the same time, in a number of ionic 
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crystals certain aspects of the phase transitions are apparently unusual from the point of view of classical physics, 
and can be explained only in a quantum-mechanical context. Such crystals contain light localized ions, and thus the 
mentioned features of their behavior point to the essential role of these light ions. In the simplest models of these 
crystals the motion of heavy complexes is ignored - their only role is to create a potential field for the light particles. 
In this case, the relevant model parameters are those which describe this field, the masses of the light particles, and the 
interaction strength. Therefore, a consistent theory of phase transitions in such models should describe the influence 
of all these parameters. We refer the reader to the survey article () for more detailed arguments in favor of the model 
given in (J2j) . 

The mathematical theory of a harmonic oscillator is quite standard, its updated and detailed presentation can 
be found in (), pages 36-44. The theory of an anharmonic oscillator with a convex potential energy is also quite 
standard. For instance, the case of H x = H x ar + bq 2r , b > and r > 2, was carefully analyzed in (). However, convex 
potentials do not have multiple wells, whereas the tunneling between the wells is the basic phenomenon in the studied 
substances. The theory outlined above in subsection III.BI among others employs quite recent results on Schrodinger 
operators obtained in (). By means of them we found out for which potentials V the domain of self-adjointness of H x 
can be established exactly. This in turn allowed us to study the gap parameter (j2Tj) in great detail. Quite often, see, 
e.g., () and (), the description of the tunneling between the wells is concentrated on semiclassical expansions for the 
gap parameter E\ — Eq, which in our case would be in negative powers of m — m p h/fi 2 - Our analysis is essentially 
different, since it covers the case of 'strong quantum effects' where semiclassical arguments are not applicable. 

2. Thermodynamic phases and the free energy 

The use of thermodynamic phases as mathematical objects is crucial in the rigorous description of thermodynamic 
properties of a given model. In classical equilibrium statistical physics, see (), (), and (), thermodynamic phases 
appear as extreme elements of the set of Gibbs measures. For quantum models described by bounded Hamiltonians, 
thermodynamic phases are pure KMS states, i.e., positive normalized linear functionals on algebras of quasi- local ob- 
servables which obey the Kubo-Martin-Schwinger conditions, see Section 6.2 in (). However, for quantum anharmonic 
crystals described by unbounded Hamiltonians, the direct construction of such KMS states appears to be impossible, 
see, e.g., the discussion in (). In this situation, the construction of thermodynamic phases as path measures realized 
in () seems to be the only possible way to settle the problem. It is performed by means of local Gibbs measures, cf. 
(|39|) . within the general scheme developed in (). The most significant point of our technique is the representation of 
Matsubara functions in the form of path integrals, see (|38|). The local Gibbs measure (f39|) which appears in this 
representation is constructed from the corresponding Hamiltonian (f2"T|) . The harmonic part X^eA ^^T* °f the latter is 
represented in (J39J) by the measure xa (|42[) , which is the Gibbs measure of noninteracting harmonic oscillators located 
in A. The interaction, the anharmonic part of as well as the external field term are taken into account in the 
energy functional (|40|) . By means of a technique developed in Chapter 2 of () the measure ([39]) is approximated by 
the Gibbs measure of a system of classical ^-dimensional oscillators living on a d + 1-dimensional lattice, see Section 
2.1 ibid. This extra dimension appears, in particular, when one approximates the integrals in (1401) by finite sums. 
Thus, the key feature of our method of studying thermodynamic states of the model © is that we describe it as a 
system of classical oscillators performing 'infinite dimensional oscillations'. 

The most significant fact about the free energy density is its universality expressed in the equality (|54")l . A deep 
fact behind it is the existence of van Hove sequences of subsets of Z d . It is known that such sequences do not exist 
for non-amenable graphs, e.g., for a Bethe lattice. In this case, the phase diagram of an Ising model is completely 
different, compared to that of this model on Z d . In particular, phase transitions can occur at nonzero h, cf. Theorem 
Q] above and (). The mentioned universality, as well as the concavity of F(h) which follows from the property ([57]) . 
allow one to get the polarization (|59[) . and hence to study phase transitions. 

3. Phase transitions, order parameters, and infrared bounds 

When one deals with thermodynamic phases, the most general definition of a phase transition is the multiplicity 
of such phases existing at the same values of the model parameters and of the temperature and an external field, 
see Chapter 7 in (). Another possibility to define a phase transition, which goes back to L. D. Landau, is to use the 
differentiability property of the free energy density F(h) as a function of the external field h. One more possibility is 
to employ an order parameter, cf. (I65|) and (I66|) . For an Ising model, all the three definitions are equivalent, see (). 
In our case, however, we managed to establish only the fact stated in Theorem[3J For versions of the model ©, phase 
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transitions were previously discussed mostly by means of the order parameter, see (). The reason for this is that, in 
such a case, one needs to control the Matsubara functions (|63| only, which can be done not only by path measures. 

The way of establishing phase transitions in quantum lattice models by means of infrared estimates was suggested 
in (). To the model ([2]), in the present form it was developed in (). The key idea of this method is to show that the 
spatial correlations do not decay to zero at large distances. This would point to the existence of a Gibbs state which 
is not ergodic with respect to the lattice translations, and hence to the multiplicity of thermodynamic phases. The 
key point here is the estimate in (|70p . where the first summand in the lower bound of D^ y is independent of x. Then 
the property in question is obtained if (|72[) is satisfied. 

B. To Section 17771 

1. Phase transitions 

As was mentioned above, in our approach the system of quantum anharmonic oscillators is described as a system 
of classical oscillators performing 'infinite dimensional oscillations'. This allows us to apply here the original version 
of the infrared estimates () adapted to the infinite dimensional case. The most transparent proof of the existence of a 
phase transition is in the |<7| 4 case. The way of passing from (f8(T[) to the key estimate (f81j) was suggested in () where 
the function / as in (|79p was introduced. The condition (|78[) is surprisingly simple. In fact, it points to a quantum 
phase transition which would occur at zero temperature. If (|T8[) holds, then there exist multiple 'ground states' of the 
model, whereas ([94]) implies that the 'ground state' is unique. Of course, such a 'ground state' ought to be defined as 
a mathematical object. For a simpler version of the model ([2]), this was done in (). The extension of Theorem [3] to 
more general models made in Theorem |4] was done by means of a comparison method elaborated in () and based on 
a tedious analysis of path measures. The result of Theorem [5] was first established in (). Its main conclusion is that a 
first order phase transition can occur without symmetry breaking. The way of getting this result was crucially based 
on the properties of path measures. 

2. Quantum effects 

As was mentioned above, the first paper where quantum effects in models like ^ were discussed is (). However, only 
the systematic use of path measures, which allows for dealing with thermodynamic states as mathematical objects, 
leads to the most complete description of such effects. Note also that a tedious analysis of the spectral properties of 
a single anharmonic oscillator was crucial. The first paper where the notion of quantum stabilization was introduced 
is (). The key parameter is the quantum rigidity ([25]) . When it is large, the oscillator 'forgets' about the details 
of the potential energy in H x in the vicinity of the origin (including instability) and oscillates as if its equilibrium 
were stable, as in the harmonic case. The results of Theorem [8] were obtained in () by means of the so called scalar 
domination inequalities obtained in (). The result of Theorem |9] was obtained in (). The first systematic study of 
phase transitions and quantum effects in the model ([2]) based on path integral methods was done in (). 
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